arXiv: 1508.0185Ivl [astro-ph.GA] 8 Aug 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 11 August 2015 (MN style file v2.2) 


Identifying mergers using non-parametric morphological 
classification at high redshifts 

Robert Thompson^’^’^, Romeel Dave^’^’^, Shuiyao Huang®, Neal Katz® 

^ NCSA, University of Illinois, Urhana-Champaign, IL 61820 
^ University of the Western Cape, Bellville, Cape Town 7535, South Africa 
^ Astronomy Department, University of Arizona, Tueson, AZ 85721, USA 
^ South African Astronomical Observatories, Observatory, Cape Town 7925, South Africa 
^ African Institute for Mathematical Sciences, Muizenberg, Cape Town 7945, South Africa 
^Astronomy Department, University of Massachusetts, Amherst, MA 01003, USA 


11 August 2015 


ABSTRACT 

We investigate the time evolution of non-parametric morphological quantities and 
their relationship to major mergers between 4 ^ ^ 2 in high-resolution cosmological 

zoom simulations of disk galaxies that implement kinetic wind feedback, i72-based 
star formation, and minimal ISM pressurisation. We show that the resulting galaxies 
broadly match basic observed physical properties of 2 ; ^ 2 objects. We measure the 
galaxies’ concentrations (C), asymmetries (A), and Gini (G) and M 20 coefficients, and 
correlate these with major merger events identified from the mass growth history. We 
find that high values of asymmetry provide the best indicator for identifying major 
mergers of > 1 : 4 mass ratio within our sample, with Gini-M 2 o merger classification 
only as effective for face-on systems and much less effective for edge-on or randomly- 
oriented galaxies. The canonical asymmetry cut of A ^ 0.35, however, is only able to 
correctly identify major mergers ^ 10% of the time, while a higher cut of A ^ 0.8 more 
efficiently picks out mergers at this epoch. We further examine the temporal correlation 
between morphological statistics and mergers, and show that for randomly-oriented 
galaxies, half the galaxies with A ^ 0.8 undergo a merger within ±0.2 Gyr, whereas 
Gini-M 2 o identification only identifies about a third correctly. The fraction improves 
further using A ^ 1.5, but about the half the mergers are missed by this stringent cut. 

Key words: galaxies:evolution, galaxies:formation, galaxies:high-redshift, galax- 

ies:structure, galaxies:general, methods:numerical 


1 INTRODUCTION 

A key to understanding galaxy evolution is the ability to 
interpret their morphologies. Galaxies are believed to grow 
and evolve via in-situ star formation, major and minor merg¬ 
ers, and gas accretion from the intergalactic and circum- 
galactic medium. Understanding a galaxy’s morphology pro¬ 
vides an important avenue for understanding how these pro¬ 
cesses are responsible for establishing the observed struc¬ 
tural features of the galaxy population across cosmic time 
(Conselice 2014), such as the origin and emergence of the 


present-day Hubble sequence. 

Visual classification along the Hubble sequence has a 
long history as a valuable descriptor of galaxy properties. 
Morphology is observed to correlate with many other prop¬ 
erties such as colour, star formation rate, surface brightness, 
and gas content. However, visual classification remains ulti¬ 
mately subjective. Hence there has been much interest in de¬ 
veloping quantitative classification schemes. Surface bright¬ 


ness profiles ( de Vaucouleurs|1948 Sersic|19M ) provide such 
an avenue for quantitatively characterising the light dis¬ 
tribution. Such parametric characterisations must nonethe¬ 
less make assumptions about the underlying functional form 
that may not hold for all stages of galaxy formation. For in¬ 
stance, they are less well-defined in the case of merging or 
other irregular galaxies. 


More recently, non-parametric methods have gained 
popularity for quantifying morphology. In this case, the 
pixelised image is directly characterised into a mathemat¬ 
ical form. Some common methods in use today are the 
Concentration-Asymmetry-Clumpiness (CAS) system (Con- 
[selice et al.|2003j , the Gini (G) coefficient, and the second- 
order moment of the brightest 20% of a galaxy’s pixels 
(M 20 ) ( |Lotz et al.|2004 ). Since these non-parametric quanti¬ 
ties do not make any assumptions about the underlying form 
of galactic structure, they are considered a less biased tool 
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to identify objects that are undergoing key transformative 
events such as mergers. 

Major mergers are thought to play an important role 
in the evolution of galaxies. In addition to transforming the 


morphology towards earlier type (e.g. Toomre & Toomre 
1972 Cox et a'r]|2008 ), they are believed to be responsible 
for generating bursts of star formation (Mihos & Hernquist 
Narayanan et al.|2010b|a ) as well as stimulating rapid 


1996 


growth of the central supermassive black hole (Sanders & 


Mirabel 1996). While it is now clear despite the resulting 


starbursts, such mergers do not dominate the cosmic star 


formation at any observable epoch (e.g. Rodighiero et al. 


2011), nonetheless their importance to galaxy transforma¬ 
tion and black hole growth makes the identification of merg¬ 
ers a key challenge in forming a complete picture for the 
emergence of the galaxy population. 

While non-parametric morphologies offer a model- 
independent approach to characterising mergers, they must 
be calibrated against other schemes, typically visual classi¬ 
fication, to interpret them. This can be most easily done at 
low redshifts, where the galaxies can be well-resolved and 
are thus straightforward to classify. In this way, CAS and 
Gini-M 2 Q have been shown to isolate mergers in particular 
regions of parameter space. Using this, CAS has been used 
to identify major merger candidates in numerous galaxy sur¬ 
veys (e.g. Abraham et al.|19'% Conselice et al. 2003| Cassata 


et al. |20C)5 Menanteau et al. 2006). Lotz et al. (2004) ar¬ 
gued that Gini-M 2 o is more effective at identifying late-stage 
mergers and classifying low-redshift galaxies. 

As galaxy surveys are extended to higher redshifts, such 
parametrisations continue to be utilised to classify mergers 


and morphologies. With surveys such as CANDELS (Grogin 
eTaTIpoTT [Koekemoer et ar]|2011| ) using Huhhle, the mor¬ 
phology of distant galaxies can be resolved at rest-optical 
wavelengths, which if interpreted properly could provide im¬ 
portant information towards understanding the evolution of 
galaxies during the peak epoch of cosmic star formation. 
However, visual morphologies of high-redshift galaxies do 
not clearly follow the Hubble sequence classification seen 
today. For example, at z ~ 2 disks tend to be thicker, more 
turbulent, and clumpier (e.g. Forster Schreiber et al.|[2Q09 


Elmegreen et al. 2009 Genzel et al. |20 11 ), ellipticals are sub¬ 
stantially more compact ( van Dokkum et al.|2014 ), and the 
fraction of irregulars increases substantially compared to to¬ 
day. It is thus not clear how well non-parametric morpholo¬ 
gies calibrated at \ow-z will apply to such distant galaxies, 
and hence whether quantitative morphology can be effec¬ 
tively utilised to characterise the high-z galaxy population. 

Cosmological hydrodynamical simulations provide a 
promising avenue to better understand how quantitative 
morphologies relate to key evolutionary events in galaxies 
such as mergers. With improving spatial resolution and a 
better understanding of the basic physical ingredients re¬ 


quired to broadly reproduce observed galaxies (Somerville 
et al.|20l5 ), they have emerged as a valuable complementary 
tool to elucidate the physical processes that shape galaxies. 
With sophisticated analysis tools, simulated galaxies can be 
processed into mock observations that may be analysed in an 
analogous fashion to real data, with the advantage that the 
true underlying distribution of the mass, light, and metals 
is known. 


Lotz et al. (2008) investigated both CAS and Gini- 


M 20 statistics for a number mock observed non-cosmological, 
isolated, merging systems in hydrodynamic simulations. It 
was found that Gini-M 2 o identified mergers at their first pas¬ 
sage and when the system finally coalesced, but did not con¬ 
sistently identify mergers at intermediate times. |Lotz et aL\ 
(2008 2010b|a ) further went on to characterise the time- 
scales for structural mergers for A and Gini-M 2 oio be of 
similar order. These simulations were used to calibrate the 
Gini-M 2 o merger identification relation for observed galax¬ 
ies out to z > 3. 


A different approach was taken by Cibinel et al. (2015), 
who classified mergers using resolved stellar maps as op¬ 
posed to single-band images. They found that a combina¬ 
tion of A and M 20 measured on the stellar mass maps of 
~ 100 HUDF galaxies provided the most accurate distinc¬ 
tion between mergers and isolated clumpy galaxies. They 
found that merger samples identified on the basis of A from 
a single H — band image resulted in a contamination as high 
as ~ 50% from clumpy galaxies, and argue that their method 
reduces this contamination. 


While informative, isolated galaxy simulations such as 
those used in these studies do not capture the cosmolog¬ 
ical context of galaxy formation. They lack crucial ingre¬ 
dients such as the frequent infall of smaller satellites (i.e. 
minor mergers), and important avenues for galaxy growth 
by smooth accretion (which is believed to be the dominant 
fuelling mode; |Keres et al.||2005 ). Given that the merger 
and accretion rates are much higher in the past than today, 
morphology could be strongly impacted by such processes. 

More recently, the focus has shifted towards using 
cosmological zoom-in simulations, in which an individual 
galaxy and its environs are evolved at much higher reso¬ 
lution relative to the entire volume. This retains the full 
cosmological context while making feasible the high reso¬ 
lution needed to study internal structural processes within 
galaxies such as violent disk instabilities ( Dekel et al.|20Q9 ). 


Snyder et al. (2015) examined the Gini-M 20 statistics for 


a number of high-resolution zoom-in galaxies at 3 > z > 1. 
They found that Gini-M 2 o was more sensitive to early stages 
of mergers {t < 0). Overall they concluded that structural 
evolution is neither universal nor monotonic, and that Gini- 
M 20 were largely sensitive to viewing angle, star formation, 
and mergers. 


In this paper we examine a number of galaxies at 
4 ^ 2 : ^ 2 within high-resolution cosmological smoothed 
particle hydrodynamic zoom-in simulations. We mock ob¬ 
serve our galaxy sample as as if they were detected in the 
H — hand CANDLES survey with HST and investigate their 
C, A, G, and M 20 values as a function of their mass histo¬ 
ries. Our primary aim is to determine the optimal approach 
for quantifying the merger history of galaxies using non- 
parametric morphological statistics. The hope is that this 
will provide a better understanding of how such statistics 
can best be used to interpret the evolution of galaxies at 
early epochs. 

This paper is organised as follows: in ^we describe our 
simulations and the methods. examines the global prop¬ 
erties of our galaxy sample. ^4]looks at the morphological 
statistics, while investigates the morphological statistics 
of each galaxy as a function of cosmic time. Finally in ^we 
summarise our conclusions. 
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2 METHODS 

2.1 Simulation Code 

Our simulations are evolved with an extended version of 


the GADGET-3 cosmological SPH code (Springel 2005 


Huang et ah 2015). It includes cooling processes using 


the primordial abundances described in Katz et ah (1996), 


with additional cooling from metal lines assuming photo¬ 
ionisation equilibrium from Wiersma et ah (2009). We ac- 


amount of artificial pressure introduced into the ISM, which 
more accurately captures the internal structure of galaxies 
while suppressing fragmentation on scales that we cannot 
resolve. 

To model star formation, we employ a molecular gas- 
based prescription. The molecular content (/H 2 ) of each 
gas particle is calculated following the analytic model of 
Krumholz et al. ( 2008[ 2009 hereafter KM T) and [McKee 
Krumholz (2010); details can be found in Thompson et al. 


count for metal enrichment from Type II supernovae (SNe), 
Type la SNe, AGB stars, and we track four elements 
(C,0,Si,Fe) individually as described in Oppenheimer 
|Dave| ( |2008| . Galactic outflows are implemented using the 
hybrid energy/momentum-driven wind (ezw) model fully de¬ 
tailed in |Dav5 et al. ( |2013| and Ford et al. (2015). 

Compared to our previous code described in |Oppen-| 


(2014). We stochastically form stars from the gas following 


a Schmidt (1959) Law as in 


p* = e*/H2—• 
tff 


(3) 


heimer & Dave (2008), we now use the ‘pressure-entropy’ 


SPH algorithm that better handles issues of two-phase in¬ 
stabilities when compared to the standard ‘density-entropy’ 
SPH formulations (see |Saitoh Makino|2013 Hopkins 2013 


for further details). Following Hopkins ( |2013 ) we use a stan¬ 
dard quintic spline kernel with 128 neighbours and adopt the 
Morris Monaghan] (1997) artificial viscosity along with a 


Balsara (1989) switch. We also implement a time-step lim¬ 
iter ( Saitoh Sz Makin^|2009 Durier Dalla Vecchia|[2012 ) 


The efficiency of star formation e* is set to 1% per local free- 
fall time tff ( Krumholz Tan|2007 Lada et al.|2010 ). This 
model allows us to regulate star formation using the local 
abundance of H 2 rather than the total gas density, which 
confines star formation to the densest peaks in the ISM as 
is typically seen in real galaxies, thereby preserving the ISM 
structure as much as the numerical resolution allows. This 
also introduces a metallicity dependence into the star for¬ 
mation efficiency, since the KMT model of molecular gas 
formation relies on the presence of dust, which may be im- 


that improves the accuracy of the time integration scheme 
in situations where there are sudden changes to a particle’s 
internal energy. 

To suppress artificial fragmentation within the inter- 


portant for suppressing early star formation (Tassis et al. 
2012 Krumholz et al.||2012 ). 


2.2 Initial Conditions 


Initial conditions are generated using the Music code (Hahn 


low their effective Jeans temperature (Katz & Gunn||1991 

& Abel|2011) assuming a flat AGDM cosmology. Gosmolog- 

Schaye & Dalla Vecchia||2008 Robertson & Kravtsov|2008 

j ical parameters are chosen to be consistent with constraints 


by ensuring that the simulation is always resolving at least 
one Jeans mass within a particle’s smoothing length. This 
translates to a threshold number density for pressurisation 
of 


3 (^kT\^ 

~ U v”g ”) V N, 


( 1 ) 


^ngb T^gas J 

where T = lO^K, uih is the mass of a hydrogen atom, k is 
the Boltzmann constant, Nngb is the number of neighbours 
within an SPH kernel, and rugas is the mass of an individual 
gas particle within the simulation (one Jeans mass here is 
NngbTUgas)- Particlcs with densities above pth and are con¬ 
sidered to be multi-phase, and a Jeans temperature floor is 
imposed 

Tjean. = (2) 

where p is the gas particle’s density. In low-resolution sim¬ 
ulations, Pth can drop below the nominal star formation 
threshold of n = 0.13 in which case we set it to 

this value; however, in the zoom simulations presented in 
this work, this does not occur. The net effect is analogous 


to adding ISM pressure (Springel & Hernquist 2003), but 


rather than directly pressurising the gas, we prevent the gas 
from cooling and thus fragmenting below the Jeans scale. 
The net effect of this implementation is in practice quite 
similar to that of Schaye Dalla Vecchia| ( |2008 ), who force 
gas particles to lie on the T oc relation, except that 
rather than using a fixed density threshold (as they use), we 
adjust our density threshold pth to naturally increase as the 
resolution of the simulation increases. This minimises the 


= 0.3, Qa = 0.7, Hq = 70, erg = 0.8, Us = 0.96. 

Our target halos are selected at z = 2 from a low- 
resolution hydrodynamic simulation consisting of 256^ dark- 
matter and gas particles in a (16/i“^Mpc)^ volume with an 
effective comoving spatial resolution of Cgrav — 1-25 h~^ kpc, 
which is All h~^ pc at z = 2. We populate each target halo 
with higher resolution particles at z = 249, with the size 
of each higher resolution region chosen to be 2.5 times the 
maximum radius of the original halo from the low-resolution 
hydrodynamic simulation. The majority of halos in our sam¬ 
ple are simulated with a single additional level of refinement 
i^grav = 625 h~^ pc or 208 h~^ pc at z = 2), while the two 
lowest mass halos are simulated with two additional levels 
of refinement {cgrav — 313 h~^ pc or 10Ah~^ pc at z = 2). 

Our six selected halos produce seven massive central 
galaxies at z = 2 that have no low-resolution particles within 
their virial radius. The galaxies span halo masses of 2x10^^ — 
lO^^M©, stellar masses of 3.6 x 10^ — 6.6 x lO^^M©, and star 
formation rates (SFRs) from 5 — 6 OM 0 /yr. We label them 
in order of increasing halo mass, and their global properties 
at z = 2 are listed in Table Note that given the stellar 
masses of our sample at z = 2 , it is likely that each of these 
galaxies will end up as an elliptical galaxy more massive 
than the Milky Way by z = 0 ( Lu et aL]|2014 2015). 


2.3 Simulation analysis with SPHGR 

We perform the simulation analysis using the recently de¬ 
veloped open-source python package Smoothed Particle 
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Table 1. Simulated Galaxy Sample oX z — 2 


Name 

Symbol 

Tfhalo 


SFR 

fgas 

G1 

■ 

1.91 X 10“ 

3.57 X 10® 

4.93 

0.54 

G2 

• 

1.98 X 10“ 

7.77 X 10® 

9.96 

0.47 

G3 

A 

3.35 X 10“ 

9.49 X 10® 

8.77 

0.57 

G4 

♦ 

4.85 X 10“ 

1.80 X 10“ 

2.51 

0.44 

G5 

0 

9.55 X 10“ 

4.05 X 10“ 

19.93 

0.39 

G6 

# 

1.15 X 10“ 

5.58 X 10“ 

37.50 

0.21 

G7 

T 

1.03 X 10“ 

6.64 X 10“ 

59.04 

0.28 


Note. — Halos used throughout this work aloug with their 
ceutral galaxies. All were selected from a low resolutiou simu- 
latiou iu a (16/i“^Mpc)^ volume. The smallest two halos (G1 
aud G2) have a comoviug spatial resolutiou of Cgrav — 
pc (208pc at z = 2), while the rest (G3-G7) have 
^grav = 625/i“^ pc (104pc at = 2). Mhaio aud M* 
are giveu iu uuits of M©, while SFRs are giveu iu uuits of 
Moyr“\ 


Hydrodyuamics Galaxy Reductiou, or SPHGR[^ ( |Thomp- 


squ|2015 ). Iu its most basic form, the code is respousible for 


ruuuiug a baryouic group fiuder to ideutify galaxies, a halo 
huder to ideutify dark matter halos, assigu galaxies to their 
respective halos, calculate halo aud galaxy global properties, 
aud fiually iterate through previous time steps to ideutify 
the most-massive progeuitors of each halo aud galaxy. lu- 
formatiou about each iudividual halo aud galaxy is collated 
iu a couveuieut, iutuitive, aud easy-to-access mauuer that 
elimiuates uuuecessary overhead iu the aualysis process. 

We use SKIeQ (Spline Kernel Interpolative Denmax) to 
identify galaxies as bound groups of baryonic particles ( | Cover- 1 
|nato et al.|1997[|Stadel|2001| ). All of the galaxies presented in this 
work are well above our nominal 64 star particle resolution limit 
( [Finlator et al.|2006| >. To identify dark matter halos and the asso¬ 
ciated baryonic particles we use t he six dimensional ph ase space 
halo finder RQCKSTAR-GALAXlEj^( [Behroozi et al.|2013| . By con¬ 
sidering particle velocities, Rockstar can much more accurately 
identify distinct groups of particles in mergers when compared to 
conventional FOF methods (e.g. [Thompson et al.|2015| >. 

After the identification of galaxies and halos via the above 
methods, SPHGR creates individual python objects for each halo 
and galaxy. These objects contain basic statistical information 
such as mass and radius, and additionally the member lists of 
gas, star, and dark matter (DM) particles (for halos) that belong 
to each object. For each SKID galaxy, SPHGR assigns halo mem¬ 
bership by cross-correlating galaxy membership information with 
halo membership to determine which halo has the most number 
of particles contributing to that galaxy. The most massive galaxy 
within each halo is classified as the ‘central’, while all others are 
classified as ‘satellites’. We calculate the Hi content of each galaxy 
using the auto-shielding approximation as detailed in|Dave et al.| 
pol^ , along with global properties such as the stellar mass, H 2 
gas mass, star formation rate (SFR), and SFR-weighted metallic- 
ity. 

SPHGR identifies the primary progenitors of each galaxy and 


^ https://bitbucket.org/rthompson/sphgr 
^ http://www-hpcc.astro.Washington.edu/tools/skid.html 
^ https://bitbucket.org/pbehroozi/rockstar-galaxies 


halo throughout previous time steps, as follows. For each halo we 
identify the halo in the previous output that contains the most 
dark matter particles from the current halo as the main progen¬ 
itor. In the case of galaxies, we utilise star particles rather than 
dark matter particles to track progenitors in a similar manner. 


2.4 Loser 

To create mock observations of our galaxies, we employ our newly 
developed code Loser (Line Of Sight Extinction by Ray-Tracing) 
to compute the emission spectrum and magnitude maps, i.e. im¬ 
ages, for each galaxy. 

Loser uses the Flexible Stellar Population Synthesis library 
(FSPS; [Gonroy et al.|[200^ [Gonroy fc Gunn||2010| > to generate 
single stellar population spectra interpolated to each simulated 
star particle’s age and metallicity. Each star’s spectrum is indi¬ 
vidually extincted based on the integrated dust column by com¬ 
puting a line-of-sight metal column density to the star towards a 
user-specified direction, using an SPH kernel-weighted integral of 
particles along the LOS. The metal column is converted to Ay 
using relations measured for the Milky Way ( |Watson]|2011| > at 
solar metallicity and above, while a metallicity-dependent dust- 
to-metal ratio is used at lower metallicities based on high-z GRB 
measurements ( jPe Gia et al. |2013t. Given Ay ^ an extinction law 
is then applied; here we use a jGalzetti et al. ] (|2000|> law. Loser 
also accounts for IGM attenuation based on [Madauj ( |1995| >, and 
nebular emission lines, though the former is unimportant for the 
bands we consider, and the latter does not play a significant role 
at the redshifts we consider. 

The attenuated stellar spectra are then grouped into pixels 
of a user-specified size. User-selected band passes are then ap¬ 
plied to obtain rest- or observed-frame magnitudes in each pixel. 
In our case, we select our pixel size to correspond to 0.03” at the 
appropriate redshift of the simulation output, which is compara¬ 
ble to the pixel scale in drizzled ( jKoekemoer et aL]|2003| > Wide 
Field Gamera 3 (WFG3) images from Hubble. The image is then 
blurred according to the diffraction limit of the incoming wave¬ 
length; e.g. for the case of F160W on WFG3 we apply a Gaussian 
smoothing to the image with a FWHM of ~ 0.05”. 

Loser accurately estimates the emergent flux in bands where 
reprocessed dust emission does not dominate, e.g. all the rest- 
ultraviolet to near-infrared bands that we consider here. Loser 
does not predict a spectrum for the far-infrared reprocessed emis¬ 
sion (though it does output the bolometric reprocessed luminos¬ 
ity) as would be predicted using a full dust radiative transfer code 
such as Sunrise ( |Jonsson||2006| >, but it is much faster and pro¬ 
vides greater user control over its input quantities such as the 
extinction law. 


3 GLOBAL PROPERTIES 

To characterise our galaxy sample and ensure that they are 
broadly in accord with observations, we first examine the evo¬ 
lution of their global properties between 4 ^ 2 : ^ 2. 


3.1 Mass histories 

We begin by examining the mass growth histories of our galaxies, 
and describe how we identify major mergers. First, galaxies are 
linked at 2 : = 2 with their most massive progenitor in previous 
outputs as described in §2.3| This allows for the extraction of each 
galaxy’s baryonic mass history, and hence we can identify at each 
output (redshift) a merger ratio IZ defined by 

( 4 ) 

M2 Mt ^ ' 
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Figure 1. Evolution of the total baryonic mass (Mgas + 
of each galaxy within our sample between redshifts 4 ^ 2 ; ^ 2. 
The end of major mergers (M 1 /M 2 ^1/4) are indicated by data 
points with black borders. 


where Mt is the mass of the galaxy at the current output, and 
galaxy’s mass in the previous output (outputs have 
an approximate spacing of At ^ 50 — 100 Myr). We consider 
merger ratios of at least 1:4 (IZ ^ 0.25) to be major mergers. We 
will not explicitly discuss minor mergers, but canonically this is 
specified by 7^ ^ 0.1. A histogram of IZ values within our sample 
separates into two components, one that peaks at 7^ ~ 0.03 and 
another that peaks at 7^ ~ 0.4, with a dip around 7Z ~ 0.25. 
Hence our major merger cut at 7^ ^ 0.25 is chosen to isolate 
this separate population of events that is distinct from just an 
extension of typical “smooth” accretion and minor mergers. 

Figure ^displays the total baryonic mass {Mgas + Af*) his¬ 
tory of each galaxy within our sample. The ending time of identi¬ 
fied major mergers are indicated by large symbols with black bor¬ 
ders. The overall mass history of these galaxies is dominated by 
relatively smooth growth ( [Murali et al.|2002[ |Keres et al.|2005| , 
but major mergers occur in all galaxies but one (G2) over this 
redshift range. There are a total of 10 major mergers that occur 
within our sample between 4 ^ 2 ; ^ 2, yielding a typical timescale 
between mergers of ~ 1.25 Gyr. 


3.2 Gas fraction 

The gas fraction {fgas = Mgas! (Af^as+Af*)) history of our galax¬ 
ies is shown in Figure [^a), again with the end of major mergers 
shown as points with black borders for reference. Overall, the gas 
fraction evolves downwards with time. As discussed in IDave et al.l 
( |2011| ), this owes to efficient outflows in small galaxies at early 
epochs that keeps galaxies from converting gas into stars. After 
major mergers, we often see a somewhat more rapid decrease in 
fgas , suggesting that these galaxies consume the gas available for 
star formation quickly during these events. 

The predicted values of fgas between 4 ^ 2 : ^ 2 for our galax¬ 
ies broadly agree with the |Erb et al.| ( |2Q06| > 2 ; = 2 observations 
shown as grey crosses, with galaxies evolving mostly down along 
the relation with time. This indicates that the ISM gas content of 
our galaxies is generally in accord with observed galaxies at these 
epochs. It is, however, somewhat higher than the gas fractions in 
the lower redshift z 1 — 1.5 PHIBBS data (shown as the grey 
contour; [Tacconi et al.|20T^ . 



Figure 2. Redshift evolution of global properties of our galaxy 
sample between 4 ^ 2 ; ^ 2. Major mergers {IZ ^1/4) are again in¬ 
dicated by data points with black borders. Panel (a) shows fgas 
for each galaxy, with observational data at 2 ; = 2 as the grey 
crosses (|Erb et al.|j2006| >, and 2 ; = I — 1.5 data from PHIBBS 
( [Tacconi et al.|2013| as the grey contour. Panel (b) shows the half¬ 
mass radius for each galaxy, the Compactness’ line for 1.4 < 2 ; < 3 
(jBarro et al.|20l^ , and 1 < 2 : < 3 observational data compiled by 
[Bruce et al.| ( |2Ql^ as the grey contour. Panel (c) shows the ratio 
of the molecular gas half-mass radius and the stellar half mass 
radius for each galaxy, the Compactness’ line for 1.4 < 2 ; < 3 
( [Barro et al.[[2013[>, and 1 < 2 ; < 3 observational data compiled 
by [Bruce et al.[ ( [2012[ > as the grey contour. Panel (d) shows the 
specific star formation rate for each galaxy with the grey con¬ 
tour representing observational data at 2 ; = 2 ( [Daddi et al.|2007| 
[Krnmholz et al.][2012[ ), and the blue contour from [Speagle et al.[ 
( [2014[ >. Panel (d) shows the specific star formation rate for each 
galaxy with the grey contour representing observational data at 


2 ; = 2 ([Daddi et al.[[2007[ 

Krnmholz et al.[[2012[), and the blue 

contour from[Speagle et al. 

(|2014||. 
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3.3 Size 

Figure Hb) shows the stellar half-mass radius (HMR*) of each 
of our galaxies between 4 ^ 2 ; ^ 2. Galaxies typically have sizes 
around 1 physical kpc at all epochs, and show no obvious trend 
with mass. The weak evolution of galaxy size with respect to 
stellar mass is in broad agreement with numerical models and 
observations out to 2 ; ~ 1 (e.g. |Somerville et al.||2008| |Brooks| 
let al.|[MTTt . Our galaxy sizes also lie within the locus of 3 > 
2 ; > 1 observations shown as the grey contour ( [Bruce et al.|201^ . 
We note that our physical resolution (i.e. gravitational softening 
length) at these epochs is generally better than ~ 200 pc, so these 
galaxies are all reasonably well-resolved. 

We also overlay the 3 > 2 ; > 1.4 line that delineates “com¬ 


pact” galaxies from [Barro et al.j ( |2013| >. Identifying the progeni¬ 
tors of massive compact galaxies has been a topic of great interest 


recently (e.g. [Barro 

et al. 

20I3[ [Patel et al.|20I3[ [Stefanon et al.[ 

MTal IWellons et al. 

[ms 

[Zolotov et al.[[20I5[ [Fang et al.||20I5[ 

Stringer et al.[[20I5 

[Graham et al.[[20I5[). In our simulations. 


galaxies grow much more slowly in half-mass radius than in stel¬ 
lar mass; in this plot they tend to move nearly horizontally with 
time. As a result, the three most massive galaxies (G7, G6, G5) 
all end up in the compact regime. G4 undergoes an early merger 
that lowers its radius into the compact regime at a low mass, but 
it quickly re-grows out of being compact, then nearly returns to 
becoming compact by 2 ; = 2. 

Mergers often greatly but temporarily increase the size. This 
may be because SKID will often identify two merging galaxies as a 
single one prior to final coalescence, resulting in a large half-mass 
extent. However, by the next output the galaxy usually settles 
back to something similar to its previous size. In detail, the fi¬ 
nal radius tends to be slightly smaller than the pre-merger one 
(though not always). Nonetheless, it is not obviously the case 
that major mergers are primarily responsible for driving galaxies 
into the compact regime. Instead, the primary driver seems to 
be simply the fact that galaxy growth in stellar mass is not ac¬ 
companied by a sizeable increase in radius at these epochs, which 
thus eventually drives galaxies into the compact regime. 

A separate question is what process removes the star-forming 
gas once it reaches the compact regime to produce red and dead 
compact ellipticals as observed ( [van Dokkum &: Brammer|2QI0] >. 
We cannot address this question directly in our simulations since 
we do not include a model for quenching star formation. Nonethe¬ 
less, if such quenching is associated with the presence of mas¬ 
sive {Mh > 10^^ Mq) halos filled with hot gas as is often spec¬ 
ulated (e.g. [Groton et al.|[^06| [Gabor &: Bournaud]|20I4| ), then 
it may be that galaxies coincidentally achieve sufficiently mas¬ 
sive halos roughly at approximately the same stellar mass where 
their radii move into the compact regime. In such a scenario, most 
quenched galaxies would start out as compact, although the exact 
fraction depends on the stochasticity in the relationship between 
quenching and halo mass (or whatever the governing parameter 
might be). We leave a more detailed study of this for future work. 


3.4 Gas radius 

Figure [^c) shows the ratio of the molecular (star-forming) gas 
half-mass radius and the stellar half-mass radius. Overall, the two 
evolve together showing, as expected, that the stars form where 
the gas is located. However, there are some interesting departures. 
Shortly after a merger, G4’s molecular half-mass radius increases 
to twice its stellar half-mass radius, before settling back down in 
size. This mostly owes to the stellar radius dropping just prior, 
suggesting a scenario where internal dynamical effects condense 
the core before gas falls back in to regrow a larger disk, which then 
forms stars allowing the stellar radius to catch up. Such effects 
have been seen in isolated merger simulations where the progen¬ 


itors are gas-rich ( [Robertson et al.|2004| as well as cosmological 
zoom simulations of early galaxy growth ( [Governato et al.|2007[ ). 

Another interesting case is G5, where just before 2 ; = 2 its 
gas radius shoots up to more than twice its stellar radius. Again, 
this happens not too long after a merger when the stellar radius 
drops a bit, suggesting a similar scenario of ongoing disk regrowth. 
This suggests that if the molecular gas radius can be measured, 
e.g. using ALMA, and compared to the stellar radius from HST^ 
outliers in this relation may identify galaxies that have recently 
undergone internal dynamical evolution via a merger, and are in 
the phase of disk re-growth. 


3.5 Specific star formation rate 

In Figure id) we plot the specific star formation rate 
(sSFR = SFR/Mi,) of our sample, with the grey contour repre¬ 
senting z 2 observational data from [Daddi et al!] ( [2007[>, and the 
blue contour representing a fit to observations by [Speagle et al.[ 

( |M4| >. 

Galaxies generally show a mild downward evolution in sSFR 
as a function of For the vast majority of the time, all our 
galaxies lie within the observed range, which evolves slowly at 
2 ; > 2 (e.g. [Krnmholz et al.|20T^ [Stark et al.|20I^ [Duncan et al.[ 

[20I4[ [Speagle et al.[[20I4[ >. At~~the final redshift, all the galaxies 

but G5 and G6 he within the range of the observations. This is a 
non-trivial success that is difficult to achieve in lower-resolution 
simulations ( [Dave[2Q0^[Sparre et al.[2QI5[[Somerville et al.[2QI5( ), 
and suggests that the higher resolution, and potentially the inclu¬ 
sion of a molecular hydrogen-based star formation prescription, 
is important for obtaining sufficiently high sSFRs at 2 : ~ 2. 

For G5 and G6, these galaxies mostly evolved within the 
observed range, but experienced fairly sudden drops in sSFR just 
prior to 2 ; = 2. This was preceded by a mild increase in sSFR 
associated with a merger event. The drop in sSFR occurs along 
with a rise in the ratio of the gas to stellar radius. Hence for 
these objects a plausible scenario is, as discussed earlier, that 
they underwent a merger and consumed their dense gas, and are 
now regrowing a gas disk that has not yet achieved a sufficient 
mass or density to restart vigorous star formation. These may 
be examples of the sort of “quenching” discussed in [Feldmann 8z\ 
[Mayer[ ( [2QI5[ >, but it is not quenching in the conventional sense 
that results in red-and-dead galaxies, since eventually the gas 
accretion will likely re-invigorate star formation (e.g. [Gabor et al.[ 

ro . - 

In summary, our galaxies lie within the range of typical ob¬ 
served galaxies at 2 : ~ 2 in terms of their gas content, size, and 
specific SFR. Hence they provide a plausible set of objects with 
which to examine morphological characteristics and their corre¬ 
lation with merger history. 


4 MORPHOLOGICAL STATISTICS 

We now examine non-parametric morphological statistics for our 
galaxies at 4 ^ 2 : ^ 2. Every galaxy in Table ^ is processed 
through Loser six times for each of our 24 outputs between 4 ^ 
2 ; ^ 2; three correspond to random orientations aligned with the 
axes of our simulation box, while the other three are composed of 
a face-on orientation, and two edge-on orientations separated by 
90°. These viewing angles will be represented by the colours cyan 
(random), blue (face-on), and red (edge-on), respectively. This 
results in 144 mock observations for each of our seven galaxies, 
or 1,008 total mock observations for our entire sample. 

We focus here on the FI60W band as the reddest WFG3 
band that provides the most direct correlation with stellar mass. 
We have also repeated this analysis using the F8I4W band of 
AGS, which falls in the rest-frame ultraviolet, but the statistics 
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and correlations we obtain are not substantially different and the 
resulting conclusions remain the same; hence we do not discuss 
F814W further here. We will discuss detailed properties of colour 
variations in future work, and confine ourselves here to the F160W 
images for calculating our non-parametric morphological statis¬ 
tics. 


4.1 Basic morphological characteristics 

To set the stage, Figure shows face-on processed images from 
Loser for each galaxy in our sample at 2; = 2. The three columns 
show the gas, star formation rate, and stellar surface densities, re¬ 
spectively, while the final two columns show the extincted Loser 
images in the WFC3 H-band (F160W) the ACS I-band (F814W) 
filters. The flux from images like these are used to calculate their 
morphological statistics in the following sections. The white bar 
at the bottom of each image represents a scale of 5 physical kpc. 

Overall, the gas surface densities clearly indicate a disk-like 
morphology, albeit with significant levels of disturbance evident 
owing to ongoing minor mergers. The SFR generally follows the 
densest gas, mainly occurring in gas with Sgas > lO^MQkpc”^ 
in accord with the |Kennicutt| |l998| > relation. The stellar compo¬ 
nent does not show as obvious a disk-like morphology, though we 
have checked that they do show rotation. The stellar map also 
shows evidence for ongoing merging activity, most strongly in the 
smaller galaxies such as G2. The gas disks are quite extended 
compared to the stars, reaching 20 kpc across. 

In the mock images, even at HS'Ps resolution, many of the 
small-scale features in these galaxies are washed out. Nonetheless, 
the fact that G2 is about to undergo a merger clearly remains 
evident; interestingly, this is the only one of our galaxies that 
does not show an identified merger from 4 < z < 2. Spiral arms 
are barely if at all evident. They are most apparent in G5, which 
resembles the grand design spiral seen by [Law et al.| ( [^12| >. 


4.2 Concentration and Asymmetry 


A basic morphological parameter of a galaxy is its light concentra¬ 
tion (C). This has long been used to morphologically characterise 
galaxies in various classification schemes (e.g. [Morgan & May-| 
|all|1957| ). It has also been used extensively to quantify trends in 
local galaxies (e.g. [Bershady et al.||2000| [Gonselice et al.||2003] >. 
Higher values of C mean that a larger amount of the galaxy’s 
light is contained within the central region. Goncentration is typ¬ 
ically defined as the ratio of light within a circular radii containing 
20% (r2o) and 80% (rgo) of the total galaxy flux ([Gonselice et al.| 
[ 200 ^ : 

C = 51og(^). (5) 

\r2o J 


The concentration is calculated using pixels within 1.5rp where rp 
is the Petrosian radius, which we define as the radius of a circular 
aperture where the flux at that radius is 20% of the flux within 
that radius (e.g., [Bershady et al.|2000| [Gonselice et al.|2003[ >. 

Asymmetry (A) is another common morphological statistic 
that quantifies the rotational symmetry of a galax y ([Abraham[ 
[et al.[[I996[ [Brinchmann et al.[[I998[ [Gonselice et al.|[2003[ ). It is 

measured by rotating the image by 180° and subtracting that 
from the original image ([Abraham et al.[[I996[ [Gonselice et al.[ 

[2000] ): 

y; 1/0 - /1801 


A = 


El/ol 


( 6 ) 


where fo is the original image pixel flux values, and /iso is the 
rotated image’s pixel flux, and the sum is over all pixels within 
1.5rp. Higher values of A correspond to higher degrees of asym¬ 
metry, which can be used to identify merger-candidates. |Gonselice| 


4>z>2 



Concentration 


Figure 4. Asymmetry A as a function of Goncentration C for our 
sample between 4 ^ 2: ^ 2. Each galaxy in our sample contains 
6 data points: 3 random viewing orientations (cyan), 2 edge-on 
views separated by 90° (red), and one face-on view (blue). Point 
size corresponds to the specific star formation rate (sSFR) of each 
galaxy, with smaller points representing smaller sSFRs. Dashed 
lines correspond to regions in which different galaxy types are typ¬ 
ically found in the local universe ([Bershady et al.|2000[ [Gonsehce[ 
let al.|2003[ . - 


let al.| l |2003||2008t , for example, set a threshold value of A ^ 0.35 
for major mergers based on their studies of nearby galaxies in 
various phases of evolution. They focused on numerous different 
data sets that included normal galaxies, starburst galaxies, dwarf 
irregulars, dwarf ellipticals, and galaxies undergoing a merger. 

C and A are traditionally plotted against each another to 
classify a galaxy’s morphology and merger state. Figure^ shows 
C vs. A for our 1,008 mock observations between 4 ^ 2: ^ 2. The 
lines indicate the region in which different galaxy types are typi¬ 
cally found in the local universe ( [Bershady et al.|2000[ [Gonselic'^ 
[et al.[2003[ ), and the point sizes are representative of each galaxy’s 
sSFR. 

Within our sample we find a median concentration of 3.3, 
and a median asymmetry of 0.49. Hence using the canonical defi¬ 
nition that identifies mergers as objects with A ^ 0.35, the median 
galaxy at 2; = 2 — 4 in our sample is a merger. Given that major 
mergers are not that frequent as shown in our mass history plots, 
this indicates that the canonical value, tuned for lower-redshift 
galaxies, is not straightforwardly applicable at higher redshifts. 
We will examine this more quantitatively later. 

Examining each of these quantities individually allows us to 
investigate potential correlations with the merger ratios extracted 
from the simulation (Q. We split our data into three groups cor¬ 
responding to the face-on, edge-on, and random viewing angles. 
Random orientation is the most directly relevant when compar¬ 
ing to the full sample of objects in a galaxy survey, but it may 
be possible to isolate a sample of face-on or edge-on objects for 
comparison to the other panels. Two power laws are fit to the 
data within each panel; one with no weighting (blue), and the 
other weighted by the galaxy’s sSFR (orange). The sSFR weight¬ 
ing can potentially assist merger identification since mergers tend 
to have higher sSFRs, and it is also information typically available 
in a multi-wavelength galaxy survey. Each panel’s legend shows 
the power-law index of each fit along with the uncertainty in the 
slope. 

Figure [^ shows the merger ratio as a function of the concen¬ 
tration of each galaxy for each viewing orientation. The horizontal 
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Log Log SsFR Log S* F160W F814W 



Figure 3. Face-on data for each galaxy from Table at z = 2. Columns 1-3 show surface density features, while columns 4 and 5 
show each galaxy processed through LOSER in the F160W and F814W bands, respectively, as described in §2.4| The white scale bar 
corresponds to 5 physical kpc. 
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dashed line shows the value of IZ = 0.25 above which we iden¬ 
tify the galaxy as having undergone a major merger since the 
immediately previous output. From this it is clear that there is 
hardly any correlation between the concentration parameter and 
our galaxy sample’s merger state. This is not entirely surprising, 
as the major merger cut criteria is usually defined in terms of A 
(Figure]^. The only significant correlation occurs in the face-on 
sSFR-weighted fit, but this shows an antz-correlation since one 
expects mergers to increase the concentration. 

Figure is constructed in a similar manner to Figure but 
for the asymmetry parameter. In contrast to the concentration, 
we find a strong correlation between asymmetry A and a galaxy’s 
merger history, with the power-law fits having positive slopes and 
fairly small errors indicating a highly significant detection of a 
correlation. However, the scatter around the best fit is still fairly 
large (~ 0.44 dex), which means that there is not necessarily a 
tight correlation between IZ and A. The correlation appears to be 
significantly strengthened by a population of high-7^ objects that 
form an upward spur at the highest A values. 

The canonical merger separation value of A = 0.35 ( |Con-| 
[selice et al.|2003t is shown as the vertical dashed line, while higher 
values of A = 0.8 and A = 1.5 are shown as dotted lines that will 
be investigated in ^ Owing to the large scatter, there are still 
many galaxies classified as major mergers with asymmetry val¬ 
ues ^ 0.35 even though their merger ratios are consistent with 
minor-mergers or in-situ growth (lower right quadrant). Out of 
all the data points within this sample, only ~ 9% with A ^ 0.35 
are classified correctly as major mergers according to their merger 
ratios extracted from the simulation. 

Two avenues can immediately be seen for potentially im¬ 
proving the fidelity of merger identification using A. First, it is 
clear that all but one of the true mergers lie at A 0.35, hence 
we can cut at a higher value of A to obtain a higher fraction 
of mergers. Second, it is possible that A identifies mergers well, 
but perhaps in an earlier phase that is not associated with the 
final coalescence, or even a later phase. We will investigate these 
avenues to improve merger characterisation using A in ^ 


4.3 Gini and M 20 

The Gini coefficient originated as an economics statistic repre¬ 
senting the income distribution of a nation’s residents. A value of 
zero means perfect equality, wherein a value of unity represents 
perfect inequality. When used in galaxy morphology determina¬ 
tions, it quantifies the “lumpiness” of the light distribution, with 
no prior assumptions about its morphology ( [Abraham et al.|2003[ 
[Lotz et al.||2Q04| >. In the context of galaxy images, a Gini coef¬ 
ficient of zero means that the flux is evenly distributed among 
all the pixels, while a value of one means that all of the flux is 
concentrated in a single pixel. 

To compute the Gini coefficient we first sort each galaxy’s 
pixels by flux in increasing order before applying the following 
equation (e.g. [Lotz et ah p004| l: 

n 

Gini=^G —TT 5^(2*- l)/i- (7) 

fn{n - 1) 

i 

Here / is the flux of each pixel, and n is the number of pixels in 
each image. 

M 20 is defined as the normalised second-order moment of the 
brightest 20% of a galaxy’s flux. This quantity is calculated by 
again sorting the pixels from highest to lowest flux, then summing 
over the brightest pixels until reaching 20% of the total flux, and 
normalising by its total second order moment (Mtot)- 

M 20 = log , while < 0.2/tot, (8) 

^ ^ 7 


4>z>2 



Figure 7. Gini — M 20 relation for our sample of galaxies plotted 
in the same fashion as Figure ^ The different regions on this 
plot represent the different morphological classifications that have 
been applied to local galaxies, and extrapolated to higher redshift 
observations. 


where ftot is the total flux of the galaxy, and fi are the fluxes in 
each pixel ordered from highest to lowest flux. Mtot is defined as 
the flux sum of each pixel multiplied by the squared distance to 
the centre of the galaxy: 

n n 

Mtot = ^ Mi = ^ fi \{xi - Xc)"^ + {yi - yc)^] , ( 9 ) 

i i 

where {xc,yc) is the galaxy’s centre, and (xi,yi) is the pixel lo¬ 
cation. 

This value encodes similar information as the concentration 
parameter ( ^4.2| >, but differs in two primary aspects: it is more 
heavily weighted by the spatial distribution of luminous regions 
because of the dependence, and it is not measured with cir¬ 
cular or elliptical apertures. [Lotz et ah] ( |2Q04| > argues that these 
differences make M20 more sensitive to merger signatures than 
the concentration parameter. 

Akin to C and A, G and M 20 are routinely plotted as a func¬ 
tion of one another which we show in Figure[^for our sample. The 
different regions on this plot represent the different morphological 
classifications that have been applied to local galaxies ( [Lotz et al.j 
|2Q08| >. Within our sample, we find a median Gini coefficient of 
0.51, and a median M 20 of -1.74. Hence our typical z = 2 galaxy 
is an Sb/Sc spiral, which is consistent with the images shown in 
Figure 

Figure 1^ plots the merger ratio (log 77.) versus the Gini co¬ 
efficient (logG), similar to Figure]^ By itself, Gini is essentially 
uncorrelated with 77., except for a weak correlation in the face-on 
orientation. Not surprisingly, the Gini coefficient alone offers a 
poor way to identify mergers. 

Figure plots the merger ratio as a function of the M 20 
parameter, similar to Figure Again, there is only at most a 
weak anti-correlation of M 20 and merger ratio, and only in the 
face-on case. Hence both Gini and M 20 in isolation only can be 
associated with mergers in the face-on case, and even then not 
reliably. 

We note that [Stott et al.j ( [20I3[ > used visual classifications 
in the HiZELS survey at 2 : ~ 0.4 to calibrate the relationship be¬ 
tween Gini, M 20 , and disturbed morphologies indicative of merg¬ 
ers. They found that a Gini-independent cut of M 20 > —1.5 was 
efficient at detecting visually-classified mergers. Our results would 
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Figure 5. The merger-ratio as a function of the concentration (C) parameter discussed in ^3.1| split into different panels corresponding 
to the different viewing angles. Point shapes correspond to different galaxies (see Table Q while point size corresponds to a galaxy’s 
sSFR. The dashed line represents the threshold for major mergers with a mass ratio of 7^ = 0.25. Two power laws are fit to each panel 
and shown with Icr error contours; the blue is a result of fitting with no weighting, and the red fit is weighted by the sSFR ( ^3.5| >. Legends 
indicate the power law index with errors. 


Face On Edge On Random 



Log Asymmetry Log Asymmetry Log Asymmetry 


Figure 6. The merger-ratio R as a function of the asymmetry (A), similar to Figure]^ The vertical dashed line represents the canonical 
asymmetry threshold of A ^ 0.35 ( [Conselice et al^|2003| , while the two dotted vertical lines represent arbitrary thresholds of A ^ 0.8 
and 1.5 ( |5.2| >. 


not support such a cut being used at 2 ; ~ 2, as it would not do 
well at identifying galaxies that have just undergone a merger. 

4.4 Merginess 

While Gini and M 20 independently do not do well at identify¬ 
ing mergers, typically merger classification relies on a combina¬ 
tion of these two parameters. To quantify how well the Gini- 
M 20 relation classifies each data point as a merger, we define the 
term ‘merginess^ as the perpendicular distance to the canonical 
dividing merger line (similar to what was done in [Snyder et al.| 


|2015| >. A point above this line (GmTM 2 o merger) is assigned a 
positive value, while a point below {Gini-M 2 o non-merger) is as¬ 
signed a negative value, given by the distance to the line. The 
numerical value of this distance in Gini-M 2 o space has no partic¬ 
ular physical interpretation, but can be used as a relative measure 
to identify how well the Gini-M2o relation classifies mergers. 

Figure shows TZ of each object in our sample as a func¬ 
tion of merginess. Every data point to the left of the vertical 
dashed line is not classified as a merger on the Gini-M 2 Q plane, 
while everything to the right is a merger. Each data point above 
the dashed horizontal line (R = —0.46) represents a true major 
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Face On Edge On Random 



Figure 8. The merger-ratio as a function of Gini, similar to Figure 


Face On 


Edge On 


Random 



Figure 9. The merger-ratio as a function of M 20 , similar to Figure 


merger (between the previous two outputs) of at least a 1:4 mass 
ratio. Points in the upper right and lower left quadrants indi¬ 
cate objects that are correctly classified; points in the other two 
quadrants are incorrectly classified. 

Figure |10| indicates that the location in the Gini- 
M 20 relation does not do particularly well at identifying mergers 
at these redshifts. Points in the upper left quadrant are major 
mergers in the simulation, yet are classified as non-mergers on 
the Gini-M 2 o plane. Points in the lower right on the other hand, 
have merger ratios consistent with minor-mergers yet are clas¬ 
sified as major mergers on the Gini-M 2 o plane. Out of the 251 
data points classified as a merger in the Gini-M 2 o relation (posi¬ 
tive ‘merginess’), only 18 (~ 7%) he above the IZ > 0.25 line and 
are a direct result of major mergers. 

Fits to the data suggest very little correlation, with the ex¬ 
ception of face-on galaxies that have a power-law slope that de¬ 
viates by more than 3cr from zero. This correlation is still much 


less significant than for the asymmetry (A). In particular, unlike 
for A where many of the true mergers lie at high-A, here the true 
mergers are not well-localised in merginess. Even for the face-on 
case, there are only 5 major mergers out of the 31 with positive 
merginess (~ 16%) that are correctly classified. 


In summary, of the non-parametric morphological statistics 
examined here, only asymmetry provides a clear correlation with 
mergers, with the Gini-M 2 Q -based merginess parameter showing 
some correlation only in the case of face-on galaxies. Recall, how¬ 
ever, that our merger quantification only picks out major mergers 
that have occurred between the previous two outputs. It is possi¬ 
ble, as was mentioned for asymmetry, that these statistics might 
not identify mergers exactly when they end, but perhaps they do 
better at identifying recent or future mergers. We examine this 
in the next section. 
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Figure 10. The merger-ratio as a function of the merginess, similar to Figure We additionally plot a vertical dotted line at 
merginess= 0 to differentiate between a merger (positive merginess) and a non-merger (negative merginess) on the Gini — M 20 plane. 


5 MORPHOLOGICAL STATISTICS IN TIME 

In this section we plot each morphological statistic and merger 
ratio over varying intervals in cosmic time. First, we qualitatively 
examine the evolution of each of our statistics, in relation to the 
evolution of the merger ratio. Then, we quantify how well each 
statistic does at tracking mergers over a range of time intervals. 
For brevity, we omit consideration of Gini and M 20 individually, 
and focus only on the merginess parameter defined by the location 
in Gini-M 2 o space. 


5.1 Evolution of statistics 

In the figures of this section, we overlay the evolution of the 
merger ratio TZ and each morphological statistic from 4 ^ 2 : ^ 2. 
We consider separately the three cases of face-on, edge-on, and 
random orientations; for the two cases with multiple data points 
(edge-on and random) we take the median value at any given red- 
shift. We plot the merger ratio 7^ as a dashed line and plot the 
dividing line between major and minor mergers as a horizontal 
dotted line; hence major mergers are when the dashed line goes 
above the dotted line. 

FigurepT]shows the concentration of each of our seven galax¬ 
ies in various orientations as a function of cosmic time. We see 
that there is no indication that G is strongly impacted by major 
mergers, in accord with our findings in §4.2| In some instances the 
concentration parameter will spike shortly before or after a major 
merger, but overall there is no consistent correlation between a 
high IZ and G. This is not surprising since G is not typically used 
to classify mergers. 

Figure ^[similarly shows the asymmetry A as a function of 
time, with IZ overlaid as before. Here we immediately notice that 
in many cases, spikes in A mimic spikes in IZ. The correlation 
is certainly not tight, but it does appear that overall when IZ is 
large, A is often large as well. Again, this refiects the results we 
obtained in §4.2| but with the additional information that there 
is a more general temporal correlation between these quantities. 

We argued that a higher cut than the canonical A ^ 0.35 
(indicated by the horizontal dot-dashed hne; |Conselice et al.|2003| 
may be more effective at identifying mergers. This is also evident 
in this plot, but the optimal cut is not obvious; each galaxy would 


require a different cut in A to allow for a complete selection of 
major mergers. For GI, G3, G4, G6, and G7 for instance, a cut 
around A > 2 would identify mergers well. But this cut would 
not work for G2 since it never encounters a major merger yet the 
asymmetry approaches two at 2 ; ~ 2.7. For G5, conversely, its 
asymmetry history stays below the [Gonselice et al.| ( |2Q03| > A ^ 
0.35 cut until the major merger event, at which point A exceeds 
0.35 and remains above the cut for ~ 3 outputs after the major 
merger. This particular galaxy has a similarly smooth growth 
history similar to G2, yet its asymmetry history is much calmer. 
G5 is the galaxy within our sample showing the most well-settled 
disk-like morphology at 2 : = 2 akin to low-redshift galaxies, which 
may help to explain why a cut of A ^ 0.35 can correctly identify 
the major merger for this galaxy alone. 

Figure |13| analogously shows the merginess as a function 
of cosmic time for the galaxies in our sample. It is clear that 
the actual merger event does not straightforwardly increase the 
merginess, as we saw in §4.4| On occasion, however, the merginess 
spikes at a time that is close to, but not exactly coincident with, 
the merger event. We will quantify this in the next section. 

5.2 A quantitative assessment of A and merginess 
in time 

We now ask the question, do A and merginess identify mergers 
that have happened within a certain time interval from the time 
of observation? In other words, if one observes a galaxy at a given 
redshift to have a high value of A or merginess, what is the prob¬ 
ability that it recently underwent or will soon complete a merger? 

To quantify this, we re-assign the values of IZ around our ten 
major merger events for TA surrounding time steps to the 7Z at 
the time of the major merger itself. This equates to adjusting the 
y-axis values in Figures and for data points ±A time steps 
around each major merger, guaranteeing that these data points 
will be considered a major merger according to their merger ra¬ 
tios. We can then measure the fraction of objects above a given 
asymmetry or merginess threshold that also have a 7Z consistent 
with a major merger. In this way, we can quantitatively answer 
the question, does a high value of A or merginess correspond to 
identifying a merger within some time interval At (corresponding 
to ±N outputs)? 
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Figure 11. The concentration parameter C as a function of cosmic time for our seven galaxies. We consolidate our 6 viewing angles 
into either face-on, edge-on, or random; because the latter two have multiple data points we take the median values. We also show the 
normalised merger-ratio history of each galaxy as the dashed line, with the dotted line representing the threshold for a major merger; a 
major merger occurs when the dashed line crosses above this dotted line. 



Figure 12. The asymmetry parameter A as a function of cosmic time for our seven galaxies, similar to Figure Asymmetry values 
above A ^ 0.35 are canonically considered major mergers ( [Conselice et al.|2003| >; we plot this threshold as a dot-dashed line. 


Figure shows these fractions for various thresholds of A 
for our different viewing orientations, as a function of the time 
interval ±At. The solid lines represent the fraction of correctly 
identified major mergers for a given At with A ^ 0.35 ( [Conseli^ 
|et al.|2003| >, while the dashed and dotted lines show A ^ 0.8 and 
A ^ 1.5, respectively. The dot-dashed line considers all values 
of A (i.e. A > 0), so effectively shows the “background noise” 
that signifies the likelihood of being within AAt of a merger at 
any given time. Error bars for the random viewing orientation 


represent the spread in the fraction of correctly identified mergers 
for the x, y, and 2 ; viewing orientations independently. 

We can see that a galaxy within our sample having A ^ 0.35 
provides little better merger identification than the noise, for any 
time interval and viewing orientation. However, increasing the cut 
to A ^ 0.8 significantly increases the ability to identify mergers. 
A galaxy satisfying this cut has a ~ 50% chance of being within 
±200 Myr of a merger. A ^ 1.5 provides even higher percent¬ 
ages, so that there is now a ~ 70% chance of a major merger 
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Figure 13. The ‘merginess’ as a function of cosmic time for our seven galaxies, similar to Figure [IT] The dot-dashed line represents the 
difference between a merger (positive merginess) and a non-merger (negative merginess) in the Gini — M 20 plane. 


within ±0.2Gyr. The percentages are somewhat higher for edge- 
on systems than face-on, indicating that A is better at identifying 
mergers in edge-on systems. 

There is also the converse question, i.e. what fraction of true 
mergers are identified as such by these cuts? For A ^ 0.35, it 
finds all but one of the true mergers (while also unfortunately 
finding many that are not). A ^ 0.8 identifies 90% (9/10) of all 
major mergers, while A ^ 1.5 identifies only 60% (6/10) of them. 
Hence, while an increased cut improves the accuracy of merger 
identification, it sacrifices the ability to identify a highly complete 
sample of merger events. 

Figure shows the analogous fractions of correctly- 
identified mergers for various merginess thresholds, as a function 
of ±At. For the random orientation case, merginess provides es¬ 
sentially no benefit for identifying mergers relative to the back¬ 
ground noise. Only in the face-on case does merginess provide 
additional information regarding a merger, and interestingly only 
for At > 200 Myr. Changing the merginess cut does little to help, 
unlike in the asymmetry case where a more stringent cut greatly 
improves the fidelity of merger identification. 

So far we have considered intervals symmetric in time, i.e. 
we have not distinguished whether a given statistic better iden¬ 
tifies pre-mergers or merger remnants. Figure shows the effec¬ 
tiveness of merger identification for A considering positive and 
negative At separately. It is clear that A preferentially identifies 
pre-coalescence mergers, and shows the same slight preference for 
identifying these in edge-on systems as before. 

Merginess, on the other hand, displays a minor preference 
for identifying merger remnants in the face-on case, as illustrated 
in Figure [T7| but edge-on and random viewing orientations still 
provide little to no benefit compared to the background noise. 
Increasing the merginess threshold typically yields a decrease in 
the fraction of correctly identified mergers when looking in the 
—At direction. In the +At direction, however, we find that in¬ 
creases in the merginess threshold provides minor benefits only 
at At > 200Myr. 

In short, a high value in the asymmetry (A), such as 0.8, 
provides the best way to identify major mergers at 2 : ~ 2 — 4. It 
fares slightly better in the edge-on systems than in face-on ones. 


but even for a random orientation a galaxy with A ^ 0.8 has a 
~ 50% chance of being involved in a merger within ±0.2 Gyr, 
and it identifies all but one merger in our sample. A higher cut 
of A > 1.5 provides a greater probability for being (temporally) 
close to a merger, at the cost of a significantly incomplete sample 
of mergers. The location in the Gini-M 2 o plane, as quantified by 
‘merginess’, is only mildly effective at identifying mergers in the 
face-on case, and ineffective in the edge-on or random orienta¬ 
tions. Moreover, no other cut in merginess is significantly more 
effective. We thus conclude that asymmetry represents a better, 
though not perfect, approach to identifying mergers in high- 2 ; 
galaxies than Gini-M 2 o . 


5.3 Resolution dependence 

To assess the resolution dependence of our results, we examine 
the lower-resolution counterparts to the galaxies listed in Table ^ 
from the original low-resolution full box simulation ( ^2.2| ). By 
subtracting the fractional value from the background noise ( ^5.2| > 
we can quantify the increase in the fraction of mergers correctly 
identified in relation to the background noise for each of our A 
and merginess cuts. 

Figures |I8| and |I9| show this fractional increase with re¬ 
spect to the background noise for the random viewing orienta¬ 
tion (we omit the face-on and edge-on viewing orientations as 
they show similar trends). The left panels show the results for 
the low-resolution galaxy sample, while the right panels shows 
the high-resolution zoom sample. As in Figures |l4|I7| the error 
bars represent the spread in the percentage increase in the frac¬ 
tion of correctly identified mergers when examining the a:, y, and 
2 ; viewing orientations independently. 

We find that the general conclusions drawn here as regards 
the efficacy of A and merginess are not overly sensitive to reso¬ 
lution. We do find that the overall fraction of correctly identified 
systems along with the background noise level both increase some¬ 
what with a decrease in resolution, but the difference between the 
two remains similar. 
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Figure 14. The fraction of major mergers correctly identified above a given asymmetry threshold when we consider the time dependent 
effects discussed in ^5.2| The solid lines represent the typical asymmetry cut {A ^ 0.35; |Consehce et al.||2003| > used to identify major 
merger candidates, while higher cuts of 0.8 (dashed) and 1.5 (dotted) are arbitrary choices to illustrate the effectiveness of A at selecting 
major mergers within our sample. The black dot-dashed line is the chance of encountering a random merger over a given as if there 
were no threshold value (background noise). Error bars represent the spread in the fraction mergers correctly identified when examining 
the X, y, and 2 ; viewing orientations independently. 



±At [Gyr] ±At [Gyr] ±At [Gyr] 

Figure 15. The fraction of major mergers correctly identified above a given merginess threshold when we consider the time dependent 
effects discussed in §5.2| The solid lines represents a merginess threshold of merginess=0 used to identify major merger candidates on 
the Gini-M 2 Q plane, while higher thresholds of 0.05 (dashed) and 0.1 (dotted) are arbitrary choices to illustrate the effect of shifting the 
merginess threshold. The black dot-dashed line is the chance of encountering a random merger over a given as if there were no 

merginess threshold (background noise). Error bars represent the spread in the fraction mergers correctly identified when examining the 
a:, y, and 2 ; viewing orientations independently. 
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Figure 16. Similar to Figure but here we plot the fraction of mergers correctly identified when separately considering the +At and 
—At time dependent effects discussed in §5.2| for the asymmetry parameter. 



At [Gyr] At [Gyr] At [Gyr] 

Figure 17. Similar to Figure but here we plot the fraction of mergers correctly identified when separately considering the +At and 
—At time dependent effects discussed in §5.2| for the merginess parameter. 


6 CONCLUSIONS 

We have performed high-resolution cosmological zoom-in simula¬ 
tions to investigate the evolution of nonparametric galaxy mor¬ 
phological statistics and their relationship to major mergers at 
2 ; ^ 2. We consider 7 galaxies with 3.5 X10® < M* < 6.6 x 10^® Mq 
and 2 x 10^^ < Mhaio < 10^^ Mq at 2: = 2 having star forma¬ 
tion rates from 5 — 60MQyr“^. By post-processing our simulated 
galaxies through Loser, we create mock images of each object as 
if it were observed with HST WFCS’s F160W filter. The fiuxes 
from our simulated images are used to calculate the concentration 


(C), asymmetry (A), Gini coefficient (G), and the normalised 
second-order moment of the brightest 20% of a galaxy’s ffux 
(M 20 ). We track the mass growth history of each galaxy from 
4 ^ 2 ; ^ 2 and identify 10 major merger events of at least a 
1:4 mass-ratio, and search for correlations with the above men¬ 
tioned morphological statistics that could most reliably identify 
such merger events. 

Our main findings are as follows: 

• The asymmetry statistic (A) strongly correlates with merger 
ratio among our simulated galaxies at these redshifts (Figures]^ 
and|12|>. Concentration, Gini, and M 20 do not substantially cor- 
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Figure 18. Fractional increase in the fraction of correctly identi¬ 
fied major mergers when compared to the background noise ( ^5.2[ ) 
for three cuts in A and random viewing orientation. Left: results 
from the low-resolution galaxy sample ( ^2.2| >, Right: results from 
the high-resolution zoom sample used throughout this work. Error 
bars represent the spread in the increase in the fraction of cor¬ 
rectly identified mergers when examining the x, y, and 2 : viewing 
orientations independently. 


Low-resolution High-resolution 

galaxy sample galaxy sample 



Figure 19. Similar to Figurefor the merginess parameter. 


relate with merger ratio. We define a ‘merginess’ parameter that 
quantifies the distance from the merger dividing line in Gini- 
M 20 space, but this also does not correlate with merger ratio, 
except for a weak correlation in the face-on case ( |4.4| ). 

• Examining A in more detail, we find that a higher thresh¬ 
old for A than the canonically-used 0.35 provides a more robust 
identification of mergers. For instance A ^ 0.8 identifies 9/10 of 
our mergers, with a much lower fraction of erroneously-identified 
mergers, ~ 50%. A ^ 1.5 further reduces erroneous identifica¬ 
tions, but only identifies about half the mergers in our sample. 

• We examine the temporal dependence of these parameters 
by quantifying how well each statistic identifies a merger that 
is about to happen or has just happened within some time At 
( ^5.2| . The merginess parameter (from Gini-M 2 o) only weakly 
identifies mergers above random sampling in the face-on case, 
but does not properly identify mergers in the random orientation 
case. (Figure]^. 

• The probability of A correctly identifying a merger increases 
with both At and the threshold value. A galaxy with A ^ 0.8 
has a ~ 50% chance of being an object that experiences a merger 
within ±200 Myr. The probability increases to ~ 70% for A ^ 1.5, 


at the cost of a less complete sample of mergers. (Figure p^. A 
also preferentially identifies pre-coalescence mergers (Figure |I^ . 

• Within our sample, compact galaxies (as defined by |Barro| 
let al.||M^ are formed naturally via in-situ growth or minor- 
mergers. In some instances major mergers can briefiy force a 
galaxy to become compact, after a period where galaxies become 
more extended as measured by their half-light radius. 

Overall, we find that the best approach for identifying merg¬ 
ers from 2 ; = 2 to 4 is to use the asymmetry statistic (A) with 
a higher cut than used for low- 2 : galaxies. With larger samples 
we will be able to quantify these statistics in more detail, and 
in future work we will also examine the sensitivity to numerical 
methodologies in greater detail. It is also possible that resolved 
colour maps may provide additional information regarding the 
evolutionary state of galaxies. It is hoped that these types of 
investigations will aid the interpretation of resolved studies of 
high-redshift galaxies that have emerged using HST and will be 
significantly furthered in the JWST era. 


ACKNOWLEDGEMENTS 

Analysis performed with pyGadgetReader ( |Thompson|2QI4| and 
SPHGR ( |Thompson|20I5| ) . Simulations were run on the Univer¬ 
sity of Arizona’s El Gato supercomputei]^ funded by NSF MRI 
grant 1228509, and the University of the Western Gape’s Timon 
and Pumbaa clusters. This work was supported by NASA grant 
NNXI2AH86G. RD and RT acknowledge support from the South 
African Research Chairs Initiative and the South African National 
Research Foundation. This publication is funded in part by the 
Gordon and Betty Moore Foundation’s Data-Driven Discovery 
Initiative through Grant GBMF456I. The authors thank J. Lotz, 
G. Snyder, and R. Somerville for helpful discussions. 

REFERENCES 

Abraham, R. G., van den Bergh, S., Glazebrook, K., Ellis, R. S., 
Santiago, B. X., Surma, P., & Griffiths, R. E. 1996, ApJS, 107, 
1 

Abraham, R. G., van den Bergh, S., & Nair, P. 2003, ApJ, 588, 
218 

Balsara, D. S. 1989, PhD thesis, , Univ. Illinois at Urbana- 
Champaign, (1989) 

Barro, G., Faber, S. M., Perez-Gonzalez, P. G., Koo, D. C., 
Williams, C. C., Kocevski, D. D., Trump, J. R., Mozena, et al. 
2013, ApJ, 765, 104 

Behroozi, P. S., Wechsler, R. H., & Wu, H.-Y. 2013, ApJ, 762, 
109 

Bershady, M. A., Jangren, A., & Conselice, C. J. 2000, AJ, 119, 
2645 

Brinchmann, J., Abraham, R., Schade, D., Tresse, L., Ellis, R. S., 
Lilly, S., Le Fevre, O., Glazebrook, K., Hammer, F., Colless, 
M., Crampton, D., & Broadhurst, T. 1998, ApJ, 499, 112 

Brooks, A. M., Solomon, A. R., Governato, F., McCleary, J., 
MacArthur, L. A., Brook, C. B. A., Jonsson, P., Quinn, T. R., 
& Wadsley, J. 2011, ApJ, 728, 51 

Bruce, V. A., Dunlop, J. S., Cirasuolo, M., McLure, R. J., Tar¬ 
get!, T. A., Bell, E. F., Croton, D. J., Dekel, A., Faber, S. M., 
Ferguson, H. C., Grogin, N. A., Kocevski, D. D., Koekemoer, 
A. M., Koo, D. C., Lai, K., Lotz, J. M., McGrath, E. J., New¬ 
man, J. A., & van der Wei, A. 2012, MNRAS, 427, 1666 

Calzetti, D., Armus, L., Bohlin, R. C., Kinney, A. L., Koornneef, 
J., & Storchi-Bergmann, T. 2000, ApJ, 533, 682 


http://elgato.arizona.edu/ 


© 0000 RAS, MNRAS 000 , 000-000 

































































































18 Thompson et al. 


Cassata, P., Cimatti, A., Franceschini, A., Daddi, E., Pignatelli, 
E., Fasano, G., Rodighiero, G., Pozzetti, L., Mignoli, M., &: 
Renzini, A. 2005, MNRAS, 357, 903 
Gibinel, A., Le Floch, E., Perret, V., Bournaud, F., Daddi, E., 
Pannella, M., Elbaz, D., Amram, P., & Due, P.-A. 2015, ApJ, 
805, 181 

Gonroy, G. & Gunn, J. E. 2010, ApJ, 712, 833 

Gonroy, G., Gunn, J. E., & White, M. 2009, ApJ, 699, 486 

Gonselice, G. J. 2014, ARA&A, 52, 291 

Gonselice, G. J., Bershady, M. A., Dickinson, M., &: Papovich, 
G. 2003, AJ, 126, 1183 

Gonselice, G. J., Bershady, M. A., & Jangren, A. 2000, ApJ, 
529, 886 

Gonselice, G. J., Rajgor, S., & Myers, R. 2008, MNRAS, 386, 
909 

Gox, T. J., Jonsson, P., Somerville, R. S., Primack, J. R., &: 

Dekel, A. 2008, MNRAS, 384, 386 
Groton, D. J., Springel, V., White, S. D. M., De Lucia, G., Frenk, 
G. S., Gao, L., Jenkins, A., Kauffmann, G., Navarro, J. F., &: 
Yoshida, N. 2006, MNRAS, 365, 11 
Daddi, E., Alexander, D. M., Dickinson, M., Gilli, R., Renzini, 
A., Elbaz, D., Gimatti, A., Ghary, R., Frayer, D., Bauer, F. E., 
Brandt, W. N., Giavalisco, M., Grogin, N. A., Huynh, M., 
Kurk, J., Mignoli, M., Morrison, G., Pope, A., &: Ravindranath, 

S. 2007, ApJ, 670, 173 

Dave, R. 2008, MNRAS, 385, 147 

Dave, R., Finlator, K., & Oppenheimer, B. D. 2011, MNRAS, 
416, 1354 

Dave, R., Katz, N., Oppenheimer, B. D., Kollmeier, J. A., &: 

Weinberg, D. H. 2013, MNRAS, 434, 2645 
De Gia, A., Ledoux, G., Savaglio, S., Schady, P., & Vreeswijk, 
P. M. 2013, A&A, 560, A88 

de Vaucouleurs, G. 1948, Annales d’Astrophysique, 11, 247 
Dekel, A., Sari, R., & Geverino, D. 2009, ApJ, 703, 785 
Duncan, K., Gonselice, G. J., Mortlock, A., Hartley, W. G., Guo, 
Y., Ferguson, H. G., Dave, R., Lu, Y., Ownsworth, J., Ashby, 

M. L. N., Dekel, A., Dickinson, M., Faber, S., Giavalisco, M., 
Grogin, N., Kocevski, D., Koekemoer, A., Somerville, R. S., &: 
White, G. E. 2014, MNRAS, 444, 2960 

Durier, F. & Dalla Vecchia, G. 2012, MNRAS, 419, 465 
Elmegreen, B. G., Elmegreen, D. M., Fernandez, M. X., & Lemo- 
nias, J. J. 2009, ApJ, 692, 12 

Erb, D. K., Steidel, G. G., Shapley, A. E., Pettini, M., Reddy, 

N. A., & Adelberger, K. L. 2006, ApJ, 646, 107 

Fang, G., Ma, Z., Kong, X., & Fan, L. 2015, ApJ, 807, 139 
Feldmann, R. & Mayer, L. 2015, MNRAS, 446, 1939 
Finlator, K., Dave, R., Papovich, G., & Hernquist, L. 2006, ApJ, 
639, 672 

Ford, A. B., Werk, J. K., Dave, R., Tumlinson, J., Bordoloi, 
R., Katz, N., Kollmeier, J. A., Oppenheimer, B. D., Peeples, 
M. S., Prochaska, J. X., & Weinberg, D. H. 2015, ArXiv e- 
prints: 1503.02084 

Forster Schreiber, N. M., Genzel, R., Bouche, N., Gresci, G., 
Davies, R., Buschkamp, P., Shapiro, K., Tacconi, L. J., et al. 

2009, ApJ, 706, 1364 

Gabor, J. M. & Bournaud, F. 2014, MNRAS, 441, 1615 
Gabor, J. M., Dave, R., Finlator, K., & Oppenheimer, B. D. 

2010, MNRAS, 407, 749 

Genzel, R., Newman, S., Jones, T., Forster Schreiber, N. M., 
Shapiro, K., Genel, S., Lilly, S. J., Renzini, et al. 2011, ApJ, 
733, 101 

Governato, F., Moore, B., Gen, R., Stadel, J., Lake, G., & Quinn, 

T. 1997, NewA, 2, 91 

Governato, F., Willman, B., Mayer, L., Brooks, A., Stinson, G., 
Valenzuela, O., Wadsley, J., & Quinn, T. 2007, MNRAS, 374, 
1479 

Graham, A. W., Dullo, B. T., & Savorgnan, G. A. D. 2015, ApJ, 
804, 32 


Grogin, N. A., Kocevski, D. D., Faber, S. M., Ferguson, H. G., 
Koekemoer, A. M., Riess, A. G., Acquaviva, V., Alexander, 
et al. 2011, ApJS, 197, 35 
Hahn, O. & Abel, T. 2011, MNRAS, 415, 2101 
Hopkins, P. F. 2013, MNRAS, 428, 2840 
Huang, S. et al. 2015, In prep. 

Jonsson, P. 2006, MNRAS, 372, 2 

Katz, N. & Gunn, J. E. 1991, ApJ, 377, 365 

Katz, N., Weinberg, D. H., & Hernquist, L. 1996, ApJS, 105, 19 

Kennicutt, Jr., R. G. 1998, ApJ, 498, 541 

Keres, D., Katz, N., Weinberg, D. H., & Dave, R. 2005, MNRAS, 
363, 2 

Koekemoer, A. M., Faber, S. M., Ferguson, H. G., Grogin, N. A., 
Kocevski, D. D., Koo, D. G., Lai, K., Lotz, J. M., et al. 2011, 
ApJS, 197, 36 

Koekemoer, A. M., Fruchter, A. S., Hook, R. N., & Hack, W. 
2003, in HST Galibration Workshop : Hubble after the Instal¬ 
lation of the AGS and the NIGMOS Gooling System, ed. S. Ar- 
ribas, A. Koekemoer, &: B. Whitmore, 337 
Krumholz, M. R., Dekel, A., & McKee, G. F. 2012, ApJ, 745, 69 
Krumholz, M. R., McKee, G. F., & Tumlinson, J. 2008, ApJ, 
689, 865 

—. 2009, ApJ, 693, 216 

Krumholz, M. R. & Tan, J. G. 2007, ApJ, 654, 304 
Lada, G. J., Lombardi, M., & Alves, J. F. 2010, ApJ, 724, 687 
Law, D. R., Shapley, A. E., Steidel, G. G., Reddy, N. A., Ghris- 
tensen, G. R., Sz Erb, D. K. 2012, Nature, 487, 338 
Lotz, J. M., Davis, M., Faber, S. M., Guhathakurta, P., Gwyn, 
S., Huang, J., Koo, D. G., Le Floc’h, E., Lin, L., Newman, J., 
Noeske, K., Papovich, G., Willmer, G. N. A., Goil, A., Gon¬ 
selice, G. J., Gooper, M., Hopkins, A. M., Metevier, A., Pri¬ 
mack, J., Rieke, G., & Weiner, B. J. 2008, ApJ, 672, 177 
Lotz, J. M., Jonsson, P., Gox, T. J., & Primack, J. R. 2010a, 
MNRAS, 404, 590 
—. 2010b, MNRAS, 404, 575 

Lotz, J. M., Primack, J., Sz Madau, P. 2004, AJ, 128, 163 
Lu, Z., Mo, H. J., Lu, Y., Katz, N., Weinberg, M. D., van den 
Bosch, F. G., & Yang, X. 2014, MNRAS, 439, 1294 
—. 2015, MNRAS, 450, 1604 
Madau, P. 1995, ApJ, 441, 18 

McKee, G. F. Sz Krumholz, M. R. 2010, ApJ, 709, 308 
Menanteau, F., Ford, H. G., Motta, V., Benitez, N., Martel, 
A. R., Blakeslee, J. P., Sz Infante, L. 2006, AJ, 131, 208 
Mihos, J. G. & Hernquist, L. 1996, ApJ, 464, 641 
Morgan, W. W. Sz Mayall, N. U. 1957, PASP, 69, 291 
Morris, J. Sz Monaghan, J. 1997, Journal of Gomputational 
Physics, 136, 41 

Murali, G., Katz, N., Hernquist, L., Weinberg, D. H., Sz Dave, 
R. 2002, ApJ, 571, 1 

Narayanan, D., Dey, A., Hayward, G. G., Gox, T. J., Bussmann, 
R. S., Brodwin, M., Jonsson, P., Hopkins, P. F., Groves, B., 
Younger, J. D., & Hernquist, L. 2010a, MNRAS, 407, 1701 
Narayanan, D., Hayward, G. G., Gox, T. J., Hernquist, L., Jon¬ 
sson, P., Younger, J. D., & Groves, B. 2010b, MNRAS, 401, 
1613 

Oppenheimer, B. D. & Dave, R. 2008, MNRAS, 387, 577 
Patel, S. G., van Dokkum, P. G., Franx, M., Quadri, R. F., 
Muzzin, A., Marchesini, D., Williams, R. J., Holden, B. P., Sz 
Stefanon, M. 2013, ApJ, 766, 15 
Planck Gollaboration, Ade, P. A. R., Aghanim, N., Armitage- 
Gaplan, G., Arnaud, M., Ashdown, M., Atrio-Barandela, F., 
Aumont, J., Baccigalupi, G., Banday, A. J., Sz et al. 2013, 
ArXiv e-prints: 1303.5076 

Robertson, B., Yoshida, N., Springel, V., & Hernquist, L. 2004, 
ApJ, 606, 32 

Robertson, B. E. Sz Kravtsov, A. V. 2008, ApJ, 680, 1083 
Rodighiero, G., Daddi, E., Baronchelli, L, Gimatti, A., Renzini, 
A., Aussel, H., Popesso, P., Lutz, et al. 2011, ApJL, 739, L40 


© 0000 RAS, MNRAS 000 , 000-000 





Identifying mergers using non-parametric morphological classification at high redshifts 19 


Saitoh, T. R. &; Makino, J. 2009, ApJL, 697, L99 
—. 2013, ApJ, 768, 44 

Sanders, D. B. & Mirabel, I. F. 1996, ARA&A, 34, 749 
Schaye, J. & Dalla Vecchia, C. 2008, MNRAS, 383, 1210 
Schmidt, M. 1959, ApJ, 129, 243 

Sersic, J. L. 1963, Boletin de la Asociacion Argentina de As- 
tronomia La Plata Argentina, 6, 41 
Snyder, G. F., Lotz, J., Moody, C., Peth, M., Freeman, P., Cev- 
erino, D., Primack, J., & Dekel, A. 2015, MNRAS, 451, 4290 
Somerville, R. S., Barden, M., Rix, H.-W., Bell, E. F., Beckwith, 
S. V. W., Borch, A., Caldwell, J. A. R., Haufiler, B., Heymans, 
C., Jahnke, K., Jogee, S., McIntosh, D. H., Meisenheimer, K., 
Peng, C. Y., Sanchez, S. F., Wisotzki, L., & Wolf, C. 2008, 
ApJ, 672, 776 

Somerville, R. S., Popping, G., Sz Trager, S. C. 2015, ArXiv 
e-prints:1503.00755 

Sparre, M., Hayward, C. C., Springel, V., Vogelsberger, M., 
Genel, S., Torrey, P., Nelson, D., Sijacki, D., Sz Hernquist, L. 
2015, MNRAS, 447, 3548 

Speagle, J. S., Steinhardt, C. L., Capak, P. L., Sz Silverman, 
J. D. 2014, ApJS, 214, 15 
Springel, V. 2005, MNRAS, 364, 1105 
Springel, V. & Hernquist, L. 2003, MNRAS, 339, 289 
Stadel, J. G. 2001, PhD thesis, UNIVERSITY OF WASHING¬ 
TON 

Stark, D. P., Schenker, M. A., Ellis, R., Robertson, B., McLure, 

R. , & Dunlop, J. 2013, ApJ, 763, 129 

Stefanon, M., Marchesini, D., Rudnick, G. H., Brammer, G. B., 
& Whitaker, K. E. 2013, ApJ, 768, 92 
Stott, J. P., Sobral, D., Smail, L, Bower, R., Best, P. N., & 
Geach, J. E. 2013, MNRAS, 430, 1158 
Stringer, M., Trujillo, I., Dalla Vecchia, C., Sz Martinez- 
Valpuesta, 1. 2015, MNRAS, 449, 2396 
Tacconi, L. J., Neri, R., Genzel, R., Combes, F., Bolatto, A., 
Cooper, M. C., Wuyts, S., Bournaud, F., et al. 2013, ApJ, 768, 
74 

Tassis, K., Gnedin, N. Y., & Kravtsov, A. V. 2012, ApJ, 745, 68 
Thompson, R. 2014, pyGadgetReader: GADGET snapshot 
reader for python. Astrophysics Source Code Library: 1411.001 
—. 2015, SPHGR: Smoothed-Particle Hydrodynamics Galaxy 
Reduction, Astrophysics Source Code Library: 1502.012 
Thompson, R., Dave, R., Sz Nagamine, K. 2015, MNRAS, 452, 
3030 

Thompson, R., Nagamine, K., Jaacks, J., Sz Choi, J.-H. 2014, 
ApJ, 780, 145 

Toomre, A. Sz Toomre, J. 1972, ApJ, 178, 623 
van Dokkum, P. G., Bezanson, R., van der Wei, A., Nelson, E. J., 
Momcheva, L, Skelton, R. E., Whitaker, K. E., Brammer, et al. 
2014, ApJ, 791, 45 

van Dokkum, P. G. Sz Brammer, G. 2010, ApJL, 718, L73 
Watson, D. 2011, A&A, 533, A16 

Wellons, S., Torrey, P., Ma, C.-P., Rodriguez-Gomez, V., Vogels¬ 
berger, M., Kriek, M., van Dokkum, P., Nelson, E., Genel, S., 
Pillepich, A., Springel, V., Sijacki, D., Snyder, G., Nelson, D., 
Sales, L., & Hernquist, L. 2015, MNRAS, 449, 361 
Wiersma, R. P. C., Schaye, J., Theuns, T., Dalla Vecchia, C., Sz 
Tornatore, L. 2009, MNRAS, 399, 574 
Zolotov, A., Dekel, A., Mandelker, N., Tweed, D., Inoue, S., 
DeGraf, C., Ceverino, D., Primack, J. R., Barro, G., Sz Faber, 

S. M. 2015, MNRAS, 450, 2327 


© 0000 RAS, MNRAS 000 , 000-000 




